B -> K(K*)X: DECAY DISTRIBUTIONS, BRANCHING RATIOS AND CP ASYMMETRIES 
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We discuss an improved theoretical description of the semi-inclusive B meson decays B — > K(K*)X. The decay distributions are 
calculated. Their branching ratios are found to be appreciable. The CP asymmetries in the neutral B modes B° — » K~ (K*~)X 
are sizable. An observation of direct CP violation and a measurement of 7 may come from these neutral B modes. 
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1 Introduction 

The semi-inclusive charmless hadronic decays B — > 
K(K*)X, where X is a hadronic recoil system containing 
no charmed and strange particles, are useful for studying 
direct CP violation and determining the weak phase 7 p. 
Direct CP violation can also be studied with exclusive or 
inclusive charmless hadronic B decays. Although theo- 
retical uncertainties in inclusive decays may be small, it 
is experimentally hard to identify final states inclusively, 
while exclusive decays have a clear experimental signa- 
ture, but the theoretical calculation is not as clean. Semi- 
inclusive hadronic decays, which lie somewhere between 
inclusive and exclusive hadronic decays, have a small the- 
oretical uncertainty and a clear experimental signature, 
thereby providing an interesting channel for studying-di- 
rect CP violation. Here[] we report on a recent studya of 
B -> K{K*)X decays. 

Charmless hadronic B — > K(K*)X decays involve 
two types of amplitudes: b — > u tree amplitudes and 
b — ► s penguin amplitudes. Direct CP violation arises 
due to the interference between two or more participating 
amplitudes with different weak and strong phases for a 
single decay mode. The effective Hamiltonian for charm- 
less hadronic B decays with AS = 1 is given by 



H eff = ^{VubKafaO! + C 2 2 + CnO n 



n—3 



+V cb V c * s J2c n O n }, 



(1) 



n—3 



where O n are local four-quark and magnetic moment 
operators and c n the corresponding Wilson coefficients. 
Both c„ and the matrix elements of O n depend on the 
renormalization scheme and scale. The effective the- 
ory based on the Wilson operator product expansion 
provides a framework to separate the short- and long- 
distance strong interaction. The short-distance strong 
interaction effects above the scale // ~ mb are incorpo- 
rated in the Wilson coefficients c n . The long-distance 
strong interaction effects below the scale [i are encoded 



in the hadronic matrix elements of the local operators 
< XK(K*)\O n \B >. 

The short-distance coefficients c„ are well known. 
They have been calculated up to the next-to-leading- 
order corrections a. Factorization has often been em- 
ployed to calculate-liadronic matrix elements. Recent 
theoretical work □BO has justified factorization in the 
heavy quark limit in the case that the ejected particle 
from the B decay is a light meson or an onium. Per- 
turbative QCD corrections to factorization can be com- 
puted in the heavy quark limit. The separation of short 
and long distance QCD using the effective Hamiltonian 
together with QCD factorization, light cone expansion 
and heavy quark effective theory leads to an improved 
theoretical description of B — > K(K*)X decays. 

2 Initial Bound State Effects 

We choose to study B° -> K~{K*-)X and B~ -> 
K°(K*°)X . The factorized matrix elements for these 
processes do not involve the transition form factors, elim- 
inating a potential theoretical uncertainty. Using the 
effective Hamiltonian (|l|) and QCD factorization □, the 
differential decay rate for B — > KX in the B rest frame 
is given by 
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2m B {2n) 3 2E K 
x / d 4 y e i y p K{d a A q/ (y)} < B\b(0)^U(0,y)b(y)\B > 



(2) 
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q = u and d for B° and B~~ , respectively. Af- 
ter factorization the long-distance QCD effects on B — * 
KX decays are contained in two matrix elements: < 
K\s^^(l — 75)9' 1 >= ifnP^ defining the decay con- 
stant f K for a K meson and < B\b(0)^U{Q, y)b(y)\B > 
incorporating the effect of the initial b quark bound 
state in the B meson. The K meson decay constant 
is known from experiment Q: fx = 159.8 ± 1.5 MeV. 
We need to know another hadronic matrix element < 
B\b(0)^U(0,y)b(y)\B >. 
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We use the light cone expansion method to handle 
the initial bound state effect. Since the dominant con- 
tribution to the decay rate comes from the light cone 
region, the light cone expansion can be used to sys- 
tematically calculate the matrix element. In the lead- 
ing twist approximation, < B|6(0)7 /3 C/(0, y)b{y)\B >= 
2P|J d£e- iiy p zl(£), where /(£) is the 6-quark distri- 
bution function Em. We then obtain the decay distribu- 
tion as a function of Ek 



dT(B -> KX) 



dE 



K 



fx E K 
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K 
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Carrying out similar calculations, we obtain the differen- 
tial decay rate for the B — ► K*X decay 



dT{B^K*X) _ 



dE 



K> 



4|«*| 2 /( 



2E K , 
m B 



(4) 



The distribution of E K ( r ) is a delta function with 
the peak at E K ( t) = mt/2 in the free quark decay 
b — > K{K*)q' in the b rest frame. Gluon bremsstrahlung 
and initial bound state effect smear the spectrum. The 
leading effect of the initial 6-quark bound state is de- 
scribed by the distribution function /(£). It is important 
to notice that the same distribution function also encodes 
the leading nonperturbative QCD effects in B ^ X^ 
u and inclusive semileptonic decays B — > Xlv Q. Uni- 
versality of the distribution function enhances predictive 
power: The distribution function can be measured in one 
process and then used to make predictions for all other 
processes. 

Several important properties Qfihe distribution func- 
tion are already known in QCD SB. When integrating 
£ from to 1, f Q d£f(£) must give 1 due to current 
conservation. If the decay can be considered to be a 
free b quark decay, then the b quark field is given by 



Kv) 



— p-iy-Pb 



6(0), one obtains 

/(0 = *(<;-— )• 

TUB 



(5) 



We can also estimate the mean < £ >= J Q d££/(£) and 

the variance a 2 = fl-d^ 2 /(£)— < £ > 2 Mfjng heavy 
quark effective theorytd. They are given bye 



mB omf y 



3m|' 



(6) 



where 



A = 



1 



4m B 



2m b 



< B\hg s G^a^h\B >, 
< B\h(iD T ) 2 h\B >, (7) 



are two parameters of heavy quark effective theory. The 
small value for a 2 implies that the distribution function 
is sharply peaked around mi, /tub- 

The detailed form of the distribution function is not 
yet known. We use the following general parametrization 
for the distribution function for our numerical analysis: 



f(0 = N 



£(i-O c 



[{£-a) 2 + b 2 } d ' 



(8) 



where N is a normalization constant which guarantees 
Jq 1 d£f(£) — 1. The four parameters (a, b, c, d) respect all 
the known properties of the distribution function. 

We have also calculated the initial bound state effect 
by directly using heavy quark effective theory. The 1 /mj, 
expansion relates the matrix element to the parameters 
of heavy quark effective theory EH. In terms of the pa- 
rameters of heavy quark effective theory, the decay rates 
are given by 



T(B - KX) = lKm b [\a\ 2 m 2 (l + 1^-^) 
»7r 6 6 
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The free quark decay b — > K(K*)q' results are repro- 
duced in the heavy quark limit. 

3 Results and Discussions 

We have calculated the decay distributions, CP-averaged 
branching ratios, and CP asymmetries in B — > K(K*)X 
decays. The phenomenological treatment has been im- 
proved with better theoretical understanding. There are 
still several sources of theoretical uncertainties, including 
the input parameters, the meson light cone distribution 
amplitudes, the distribution function, the renormaliza- 
tion scale, and the unknown power corrections. All the 
results obtained should be understood as valid within the 
theoretical uncertainties. 

The distribution of the K(K*) energy calculated in 
the heavy quark effective theory approach is the same 
delta function with the peak at E K {,) = m b /2 as in 
the free quark decay b — ► K(K*)q'. The light cone ex- 
pansion approach is capable of accounting for the initial 
bound state effect on the decay distribution. The re- 
sulting spectrum spreads over the full kinematic range 
< E K {,) < ms/2 and depends strongly on the form of 
the distribution function. However, if the cut E K {,) > 2.1 
GeV is applied to suppress the background, more than 
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Figure 1: Kaon energy spectrum in B° — t K~ X. The solid and 
dashed curves correspond to two very different forms of the distri- 
bution function /(£), respectively. 



0.84) x 1(T 4 [(0.67 - 0.74) x 10" 4 ] for B~ -> K°(K*°)X, 
depending on the value of the CP violating phase 7. In 
the heavy quark effective theory approach, the branching 
ratios are decreased by about 10% and the CP asymme- 
tries are not affected in comparison with the free 6-quark 
decay. The three estimates (free quark decay approxi- 
mation, light cone expansion and heavy quark effective 
theory method) all give the same order of magnitudes for 
the branching ratios and CP asymmetries. The branch- 
ing ratios for B — > K(K*)X are of order 10~ 4 and the CP 
asymmetries in the neutral B modes B° — > K~ (K*~)X 
are sizable and can be measured at the B factories. 
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97% of events will survive. For illustration, the kaon 
energy spectrum in B° — > K~X is shown in Fig. 1, 
computed in the light cone expansion approach assuming 
7 = 60°. 

One can also calculate the hadronic invariant mass 
spectrum dT/dMx in B — > KX . Mx is related with Ek 
in the B rest frame through M x = (Pg — Pk) 2 = rn B ~ 
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Figure 2: CP-averaged branching ratio for B° — > X. 



Figure 6: CP asymmetry in B° — » K X. 





Figure 3: CP-averaged branching ratio for B 



K°X. 



Figure 7: CP asymmetry in B — ► K°X. 
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Figure 4: CP-averaged branching ratio for B° — > if* X. 



Figure 8: CP asymmetry in B° -> K*~X. 
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Figure 5: CP-averaged branching ratio for B 



K*°X. 



Figure 9: CP asymmetry in B 



K*°X. 
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